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Abstract 

Three-body resonances in the KNN system have been studied within a framework of the 
KNN — irYN coupled-channel Faddeev equation. By solving the three-body equation, the en- 
ergy dependence of the resonant KN amplitude is fully taken into account. The S-matrix pole has 
been investigated from the eigenvalue of the kernel with the analytic continuation of the scattering 
amplitude on the unphysical Riemann sheet. The KN interaction is constructed from the leading 
order term of the chiral Lagrangian using relativistic kinematics. The A(1405) resonance is dy- 
namically generated in this model, where the KN interaction parameters are fitted to the data of 
scattering length. As a result we find a three-body resonance of the strange dibaryon system with 
binding energy, B ~ 79 MeV, and width, T ~ 74 MeV. The energy of the three-body resonance is 
found to be sensitive to the model of the 1 = KN interaction. 
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I. INTRODUCTION 

The analysis of the kaonic atoms [l| revealed an attractive A'-nucleus interaction. Al- 
though the strength of the attraction depends on the parametrization of the density de- 
pendence of the optical potential [1] and the theoretical study of the K optical potential 
suggests a rather shallow potential 2J, there has been a great interest in the possibilities 
of .ft'- nucleus bound states in recent years. Akaishi and Yamazaki j^, 4] studied the kaon 
bound states in light nuclei and found deeply bound kaonic states, for example, B ~ 100 
MeV for |.H. In their study, the kaonic nuclear states were investigated by using the K 
optical potential, which is constructed by folding the g matrix with a trial nuclear density. 
The potential model of KN — 7rS interaction is determined to reproduce the A(1405) and 
the scattering length. The kaonic nuclear states are further studied by using a method of 
antisymmetrised molecular dynamics [s| using the KN g matrix. 

Among the simplest K nucleus state, the K~pp state, which has strangeness S = — 1, total 
angular momentum and parity J 71 " = 0~, and isospin 1 = 1/2 dibaryon state, is expected to 
have largest component of the 1 = KN. An experimental signal of the K~pp bound state is 
reported by the FINUDA Collaboration from the analysis of the invariant mass distribution 
of A — p in the K~ absorption reaction on nuclei [6|. The reported central value of the 
binding energy, B, and the width, T, are (B,T) = (115,67) MeV, which is below the nUN 
threshold energy. This data may be compared with the predicted values (B, T) = (48, 61) 
MeV in Ref. |4J. However it was pointed out that the data can be understood by the two- 
nucleon absorption of K~ in nuclei together with the final state interaction of the outgoing 
baryons Q]. 

In the attractive interaction of kaon in nuclei, the resonance A (1405) in the s-wave and 
1 = channel KN scattering state plays an essential role. The energy of the A(1405) is 
below the KN threshold and strongly couples with the 7rS state. Although the kaonic nuclear 
states have been studied so far by using the KN g matrix or optical potential, it might be 
very important to examine the full dynamical calculation of KN — 7rE system by taking into 
account the energy dependence of the resonance t matrix and the coupling with the KN — tiH 
channel explicitly. Such a theoretical study may be possible in the simplest kaonic nuclei 
with baryon number B = 2 system. In this work, we study the strange dibaryon system by 
taking into account the three-body dynamics using the KNN — nHN — nAN (KNN — nYN) 
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coupled-channel Faddeev equation with relativistic and non-relativistic kinematics. 

Methods to investigate resonances in the three-body system have been developed in the 
studies of the three-neutron Q, [9} , irNN dibaryon [h3, [h land YiNN hypernuclei 1(3, ]j[ Q . 
In this work, we employ a method started by Glockle 8j and Moller |9j and developed by 
Matsuyama and Yazaki, Afnan, Pearce and Gibson 1(3, 
in the unphysical energy plane from the eigenvalue of the kernel of the Faddeev equation. 
To analytically continue the scattering amplitude into the unphysical sheet, the path of 
the momentum integral must be carefully deformed in the complex plane to avoid possible 
singularities. 

The most important interaction for the study of the strange dibaryon system is for the 



KN states. The internal 



structure of the A(1405) has been a long standing is- 
sue. The chiral Lagrangian approach can describe well the low energy KN 
reaction with the meson-baryon dynamics. A genuine q 3 picture of the A (1405) coupled 
with meson-baryon [l7] may not yet be excluded. Though previous studies of the KNN 
system used phenomenological models of the KN potentials, we use s-wave meson-baryon 
coupled-channel potentials guided by the lowest order chiral Lagrangian. With this model, 
the strength of the potentials and the relative strength of the potentials among various 
meson-baryon channels are not parameters but are determined from the SU(3) structure of 
the chiral Lagrangian. In this model, the A(1405) is an 'unstable bound state', whose pole 
on the unphysical sheet will become the bound state of KN when the coupling between the 
KN and the irY is turned off. We examine a relativistic model as well as a nonrelativistic 
model to account for the relativistic energy of pion in the ttYN state. 

We briefly explain our KNN — ttYN coupled-channel equations and the procedure to 
search for the three-body resonance in section 2. The model of the two-body interactions 
used in this work is explained in section 3. We then report our results on the KNN 
dibaryon resonance in section 4. This work is the extension of the early version of our analysis 
reported in Ref. [18 ]. Recently Shevchenko et al. 19] performed a similar study of the KNN 
system using Faddeev equation starting form the phenomenological KN interaction within 
a nonrelativistic framework. The comparison of our results with theirs will be discussed in 
section 4. 
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II. COUPLED CHANNEL FADDEEV EQUATION AND RESONANCE POLE 



We start from the Alt-Grassberger-Sandhas(AGS) equation 2(| for the three-body scat 



tering problem. The operators £/y of the three-body scattering satisfy the AGS equation 

U itj = (1 - 5 itj )G 1 + y]t n G U nt j. (1) 

Here we label the pair of particles j, k by the spectator particle i = 1,2,3. The two-body 
t matrix tj of particles j, k with the spectator particle i is given by the solution of the 
Lippmann-Schwinger equation: 

U = v i + ViG ti. (2) 

Here Go = 1/(W — H + ie) is the free Green's function of the three particles, and W is the 
total energy of the three-body system. 

When the two-body interactions i>, are given in separable form with the vertex form factor 
\gi > and the coupling constant ji as 

Vi = \9i >ii< g%\, (3) 
the AGS-equation of Eq. (JT]) is written in the form 

X it j(pi,pj,W) = (1 - S itj )Z itj (p i ,p j , W) + J dp n Zi tn (p i ,p n ,W)T n (W)X nd (p n ,p j ,W). 



(4) 



The amplitude X it j is defined by the matrix element of Ui j between state vectors Go\pi, gi > 
as 

Xi,j(Pi,Pj,W) = < p^gilGoUijGolp^gj > . (5) 

The state vector \p i: gi > represents a plane wave state of the spectator i and the state vector 
\gi > of the interacting pair. 

The driving term Z it j of Eq. (j4]) shown in Fig. QJa) is given by the particle exchange 
mechanism defined as 



Zi,j(pi,Pj,W) = <p h gi\G G \pj,gj> (6) 

(7) 
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FIG. 1: Graphical representation of (a) one particle exchange interaction Zij(pi,pj,W) and (b) 
two-body t-matrix Tj(W). The relative momentum of the interacting particles is given by % for 
spectator particle i. 

Here the momentum of the exchanged particle i,j) is given as p k = —pi — Pj and g((ji) 
is the vertex form factor of the two-body interaction g{q,j) =< gi\qi >. The energy Ei(pi) is 
given by E^pi) = rrii +pf/2rrii for the nonrelativistic model and E^pi) = >/ntf + pf for the 
relativistic model. The relative momentum is given by ^ = (rrikPj — rrijP^) / {rrij + m*,) for 
the nonrelativistic model, while we define = |<fj| for the relativistic model as 



,W? + Tn3 — in 
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2W t 



3 k\9 2 

r — ml 



(8) 



wi = ^(EjW+Eumr-fr (9) 

The two-body t matrix can be solved for the separable interaction as 

U = \gi > Ti(W) < gt\. (10) 
Here the 'isobar' propagator r iy illustrated in Fig. [T](b), is given as 

n{W) = [i/ 7i - / d% - Idii^Qi) P i-i 

7 PF - ^(pi) - E jk (pi,qi) 

The two-body t matrix depends on the energy E^pi) of the spectator particle. Here Ej^ is 
the energy of the interacting pair given as Ejk(pi, qi) = rrij + rrik + pf / (rrij + rrik) + qf / ' \i% for 



the non-relativistic model, while Ejk(pi,qi) = a/ (Ej(qi) + Ek(qi)) 2 + pf for the relativistic 
model. The reduced mass is defined as fa = rajm^^rrij + rrik). 
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Following the standard method of angular momentum expansion 21(, the AGS equation 
reduces to following coupled integral equations by keeping only s-wave states: 



n J 

x K itn {p i ,p n ,W)X n j{p n ,p j ,W). 



(12) 



Here we used a simplified notation for the kernel K = Zt, which can be written as 

g*(qi)g(qn) 



Ki,n(Pi,Pn, W) = 2vr / d(pi ■ p 



r n (W). (13) 



W - Ei(pi) - Ej{p n ) - E k (pi + p n ) 

The formulas given above are valid for the spinless and distinguishable particles without 
channel coupling among the Fock-space vectors. In our KNN resonance problem, we have 
included the following KNN and ttYN states: 



\a> = \Ni,N 2 ,K 3 >, 
\b> = \Ni,Y 2 ,ir 3 >, 
\c> = \Y h N 2 ,7r 3 >, 



(14) 
(15) 
(16) 



with Yj is Ej or Aj. After anti-symmetrizing the amplitude for identical particles of nucle- 
ons [ll|, we obtain the following forms of the coupled AGS equations, 



J Y K ,Y K 



Zd.Y K 
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X, 



d,Y K 



(17) 



V X N*,Y K J 



Here we have suppressed the spin-isospin quantum numbers, the spectator momentum pj 
and the total energy of the three-body system W in Z, X and r for simplicity. The concise 
notation of Yk, Y^, d and N* represents the 'isobars' and their decay channels. The decay 
channels of isobars Y K , Y^, d and N* are KN(I = 0,1), 7rE(I = 0,1) and nA(I = 1), 
NN(I = 1) and ttN(I = 1/2,3/2), respectively. Here I is isobar isospin. Those indices 
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uniquely specify the three-body states of X and Z except N* showing EiV* and AiV*. 
Therefore we have a nine-channel coupled equation of Eq. (ITT)) for spin singlet, s-wave three- 
body system. The explicit form of Eq. (JT)) when we include spin-isospin is summarized in 
the Appendix. 

The dominant Fock space component is expected to be \KNN >, and therefore the 
most important amplitudes are Xy Kj y K and Xd,y K - They couple to each other through 
the kaon exchange Zy K y K and nucleon exchange Zy K ,d mechanisms. Notice, however, the 
■kYN component is also implicitly included in Ty K y K when we solve the two-body KN — 
ixY coupled-channel equations. The tcYN components, Xy W) y K and Xn*,y k , couple with 
the KNN components through the pion exchange mechanism Z^*^ and the ttN and ttY 
'isobars' Tn*,n*i Ty nt y. The pion exchange mechanism may play an important role in the 
width of the resonance. In this work, we have not included weak YN interaction. It was 
found in Ref. 19| that the YN interaction plays rather minor role in this strange dibaryon 



system. 

To find the resonance energy of the three; 
TTJ), we follow the method used in Refs 



three-b ody system 

. fly, 000. 



j pdy system using the AGS equation of Eq. 

The AGS equation of Eq. ( TT3T) 



is a Fredholm-type integral equation with the kernel K = Zr. Using the eigenvalue r] a (W) 
and the eigenfunction |0 a (W / ) > of the kernel for given energy W, 

Zr\cj> a {W) > = r} a {W)\(j) a {W) >, (18) 

the scattering amplitude X can be written as 

" \ T^jw) ' (19) 

At the energy W = W p where f] a {W p ) = 1, the amplitude has a pole, and therefore W p gives 
the bound state or resonance energy. 

Since a resonance pole appears on the unphysical energy Riemann sheet, we need analytic 
continuation of the scattering amplitude. We use here the nonrelativistic model to explain 
a method of analytic continuation, which is based on Refs. 0, 0] . At first we examine the 
singularities of the kernel of Eq. (1131) . Above the threshold energy of the three-body break 
up W > mi + rrij + rrik, Z(pi,p n , W) has logarithmic singularities. The branch points appear 
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Imp n 



Pn 



FIG. 2: The singularities of the one particle exchange interaction Z(pi,p n ,W) in the complex p n 
plane at W = E + ie and the real pt. 



at p n = ±p Zlt2 , where 



^' 2 

i 

Pz., = Jll'Vn,-^, 

Vj 



p Zl = -^Vi + x p>HWtn--vt 



(20) 
(21) 



with 



m(m k 
rrii + mk 

mi + rrij + m k ' 
W 



rrii 



For given > 0, the cuts run from to pz 2 above the positive real axis of complex p n 
plane and from — pz x to —pz 2 below the negative real axis as shown in Fig. [2j while the 
integration of momentum p n in Eq. ( JT3l) is along the real positive axis. 

Let us consider the case when W has a negative imaginary part. For given > 0, the 
cut from pz x to pz 2 moves into the fourth quadrant across the integration contour of p n . 
Assuming the integrand of Eq. f|T2|) is an analytic function around real positive p n , one can 
perform an analytic continuation of the amplitudes by deforming the integration contour 
along the logarithmic singularity as shown in Fig. [3] and then we obtain amplitudes on the 
unphysical Riemann sheet. 
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FIG. 3: The integration contour C and the singularity of Z at TV = E — iT/2 and real value of pi. 

In principle it might be possible to solve the AGS equation keeping the momentum vari- 
ables real and taking into account the discontinuity across the cut. The moving logarithmic 
singularities depending on pi make it difficult to solve the integral equation. To overcome 
this problem we deform the integration contour of pi, p n , into the fourth quadrant of the 
complex momentum plane so that we take into account the contribution of the cuts. As an 
example of our KNN — ttYN problem, we choose the integration contour of p n as shown in 
solid line in Fig. HI Here we take the energy W = 10 — z35 + + + MeV, which 
is below the mass of KNN and above the irYN. The shaded region in Fig. H] shows the 
cuts of 'Z' for the pion exchange mechanism. The cuts become 'forbidden regions' because 
the position of the cuts depends on pi, which runs the same integration contour as p n . In 
our numerical calculation, we studied all the 'forbidden regions' for n, N and K exchange 
mechanism and determined the integration contour. With the integration contour C in Fig. 
IU we choose the physical sheet of KNN. 

The singularities of the isobar propagator t(W) arises from the three-body Green's func- 



tion in the integrand of r. The poles are at q n = ±^J2^jW t h — jpPi- Since q s = (pz 1 +Pz 2 )/2, 

we can analytically continue it into the same unphysical sheet as the case in Z as long as 
we keep the same deformed contour as the one used in Z. Another singularity we have to 
worry about is the singularity due to the two-body resonance. Since our KN — 7rS system 
has the two-body resonance A(1405), the cut starts from the two-body resonance energy in 
the complex energy plane. To examine this, we write the approximate energy dependence 
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FIG. 4: The logarithmic singularities of the tt exchange mechanism Z(pi,p n , Z) at Wth = 10 — «35 
MeV in the complex p n plane. C is integration contour of p n and pi. 



of the r as 



n{W) 



(22) 



Here pn and are the momentum and mass of the spectator nucleon. The reduced mass 
of the spectator nucleon with the isobar pair KN or 7rS is denoted as t]n and E A * is the 

2 

pole energy of A (1405). At W = + E\* + with p^ on the contour C in Figs. [5] (b), 
the two-body t-matrix has a singularity, which is plotted as a solid line in Fig. 0(a). We 
illustrate the typical trajectories of the three-body resonance pole W = W p as curves A and 
B in Fig. If the pole trajectories A and B intercept the two-body A r A(1405) cut, then 
the analytic continuation to the A^A(1405) unphysical energy sheet must be examined. The 
same situation from the Pn plane is shown in Fig. [5] (b). The momentum p* corresponding 
to the energy W p of the three-body resonance is determined by 



V 



±J2r] N (W p - E A , - m N ) 



(23) 



If p* intercepts the contour C, we have to take care of the analytic continuation of the 
A^A(1405) energy sheet. As will be seen in section 4, the trajectories of the three-body 
resonance in our calculation follow line A of Fig. 0(a) and do not intercept the singularity 
of the two-body resonance. 
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FIG. 5: The singularities due to the three-body resonance and the A(1405) in (a) the complex 
energy plane and (b) the momentum plane. 

III. MODEL OF THE TWO-BODY INTERACTIONS 

We take into account the KN interactions in J n = l/2~, 1 = and 1 = 1 states, the 
irN interactions in J 71 = l/2~, 1=1/2 and 3/2 states and the AW interaction in I = 0, 1 So 
state. Our s-wave meson-baryon interaction is guided by the leading order effective chiral 
Lagrangian for the octet baryon ip B and the pseudoscalar meson fields given as 



Lint — 



8F3 



tr(^7 M [[0,^0],^]). 



(24) 



The meson-baryon potential derived from the chiral Lagrangian can be written as 

-> o\ir i - ^ 1 m + m a 

< V , p\Vbm\P, a > = ~Ci3 >a 



(25) 



Here p and ff are the momentum of the meson in the initial state a and the final state 
(5. The strength of the potential at zero momentum is not an arbitrary constant but is 
determined by the pion decay constant F n . The relative strength between the meson-baryon 
states is controlled by the constants Cp >a which are basically determined by the SU(3) flavor 
structure of the chiral Lagrangian. The parameter of our model is the cutoff A of the 
phenomenologically introduced vertex function g a {p) = A^/^p 2 + A^) 2 . 
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The most important interaction for the study of the KNN system is the 1 = KN 
interaction. We describe the KN interaction by the coupled-channel model of the KN and 
the 7rS states. The constants Cp >a for this channel are given as Cj^n-kn — ^i^kn—kX = 
— and C^-ttT, — 8. The cutoff A is determined by fitting the scattering length a 1 ^ = 
— 1.70 + i0.68 fm of Ref. 22j]. The values of A are around 1 GeV and are given as model 
(a) in tables [T]and [II] for the nonrelativistic and the relativistic models. In general, the form 
factors of the relativistic models are hard compared with those of the non-relativistic models 
because of the weak relativistic kinetic energy. We found a resonance pole at W = 1420 — i30 
MeV for the non-relativistic and the relativistic models. The relativistic kinematics might 
be important in describing nY channel because of the small pion mass. We choose this 
model ( standard parameter of the KN interaction. 

The KN scattering lengths are not very well constrained from the data. The ranges of the 
KN scattering lengths are studied within the chiral unitary model in Ref. 23) . In this work, 
we simply examined models with the scattering length a 1 ^ = (— 1.70±0.10) +i(0.68±0.10) 
fm in order to examine the sensitivity of the energy of the three-body resonance on the input 
model of the two-body interaction. The cutoff A's for those models are given as models (b)- 
(e) of Tables HI and HT1 The values of the resonance energy are about 1415 ~ 1425 MeV, and 
the width 50 ~ 70 MeV, which are close to the values of the chiral model in Ref. [if]]. One 
can notice that there is a correlation between the real (imaginary) part of the pole energy 
of the A (1405) and the imaginary (real) part of the scattering length. Those resonance 
energies are slightly larger than the pole energy reported in Ref. 24j. Therefore as a last 
model, model (f) reproduces the deeper resonance energy 1406 — i25 MeV of Ref. 24j. The 
scattering length of this model is —1.72 + i0.44 fm, which is, however, somewhat different 



from the value — 1.54 + z0.74 fm in Ref. 



24]. 



The 1 = 1 KN interaction is described by the KN — 7rS — 7rA coupled- channel model. 



The coupling constants Cp :0l are C 



KN-KN 



2,C 



KN- 



KN-ttA 



-V6, CVe-ttE = 4 



and CVe-ttA = CVa-ttA = 0. The cutoff A's are determined to fit the imaginary part of the 
scattering length of Ref. 22] . which are given as model (A) in Tables [TTT1 and HVl for the non- 



relativistic and the relativistic models. The rea' 

0.37 + i0.60 fmof Ref. 



part of the scattering length of those models 



is larger than a 1 ^ 
model (aA) , which is model (a) for 1 = and model (A) 



221 ] . The K p scattering length predicted from 



the central values of the two kaonic hydrogen data 25 



26 



or / = 1 interactions, is between 



271 ] . To study the sensitivity of 
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the models of I = 1 KN interaction to the resonance energy of K~pp system, we constructed 
model (B) given in Tables III II and IIV1 A similar model of KN interaction is developed to 
study K~d scattering 28]. The range of the vertex form factor found in in Ref. [28], which 
is monopole form factor with 880MeV cutoff mass, is comparable to ours. 

The total cross sections of K~p reactions predicted from our models (aA), (aB) and (fA) 
are shown in Fig. [6] together with the data 



29 



30 



31 



32 



331 ] . The models (aA) and (aB) 
describe well the K~p — > K~p(Fig. [6h), K~p — > 7r + E~(Fig. Eb) and K~p — > 7r~£ + (Fig. 
Efc) reactions, where both 1 = and 1 = 1 interactions contribute to the cross section. The 
models of / = (I = 1) can be tested from K~p 7r°S°(Fig. EH) (K~p 7r°A(Fig. 
Efe)) reactions, where models 'a' and 'A/B' describe the cross sections well. The model (fA) 
tends to give smaller cross sections. It is noticed, however, as we will see, that the resonance 
energy of the K~pp system is more sensitive to the 1 = KN interaction and less sensitive 
to I = 1 interactions, while both 1 = and 1 = 1 interactions are equally important to 
describe K~p cross sections and kaonic hydrogen data. 





KN(MeV) 


TrS(MeV) 


Scattering Length (fm) 


Resonance energy ( Me V) 


(a) 


1095 


1450 


-1.70 + ^0.68 


1419.8 - i29.4 


(b) 


1105 


1550 


-1.60 + ^0.68 


1422.2 - i33.7 


(c) 


1085 


1350 


-1.80 + ^0.68 


1418.5 - i25.0 


(d) 


1120 


1340 


-1.70 + ^0.59 


1414.6 - z29.4 


(e) 


1070 


1540 


-1.70 + *0.78 


1424.3 - i28.3 


(f) 


1160 


1100 


-1.72 + *0.44 


1405.8 - i25.2 



TABLE I: The cutoff parameters, scattering length and the resonance pole of the relativistic models 
of 1=0 KN — ttT, interaction. 



The form of the s-wave irN interactions is taken as Eq. (f25|) . The constant C a ^ is 4 
for 1 = 1/2 and —2 for 1 = 3/2 states. The parameters of the potentials are determined 
by fitting the scattering length and the low energy phase shifts. For 7=1/2 state, the 
strength of the potential is modified as \Cp, a by introducing a phenomenological factor A to 
describe the data of the scattering length (0. 1788 + 0. 00h0)m~ 1 [34j and the phase shifts [35]. 
The fitted parameters A and A are shown in Table |V] together with the scattering length 
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KN(MeV) 


vrS(MeV) 


Scattering Length (fm) 


Resonance energy ( Me V) 


(a) 


946 


988 


-1.70 + *0.68 


1420.1 - i30.1 


(b) 


954 


1035 


-1.60 + ^0.68 


1422.4 - i34.7 


(c) 


940 


944 


-1.80 + *0.68 


1418.7 -i26.0 


(d) 


968 


933 


-1.70 + ^0.58 


1414.3 - i30.5 


(e) 


927 


1031 


-1.70 + *0.78 


1424.7 - i29.0 


(f) 


1000 


800 


-1.72 + ^0.43 


1404.8 - i25.5 



TABLE II: The cutoff parameters, scattering length and the resonance pole of the nonrelativistic 
models of 1=0 KN — ttT, interaction. 





KN(MeV) 


vrS(MeV) 


vrA(MeV) 


Scattering Length(fm) 


(A) 


1100 


850 


1250 


0.68 + iO.60 


(B) 


950 


800 


1250 


0.65 + i0.46 



TABLE III: The cutoff parameters, scattering length of the relativistic models of 1=1 KN — ttY 
interaction. 



calculated using the models. The models describe well the Su phase shifts up to 1.2 GeV 
as shown in Fig. [71 

For the 1 = 3/2 ttN scattering, the irN potential is constructed so as to reproduce the 
scattering length (—0.0927 ± O.OOOS)??!^ 1 [34| and the S31 partial wave phase shifts data. 
Here we introduced a modified dipole form factor as 

A 4 



g(p) 



(f + A 2 ) 5 



x (1 +af). 



(26) 





KN(MeV) 


vrS(MeV) 


vrA(MeV) 


Scattering Length(fm) 


(A) 


920 


960 


640 


0.72 + f0.59 


(B) 


800 


940 


660 


0.68 + i0.45 



TABLE IV: The cutoff parameters, scattering length of the nonrelativistic models of 1=1 KN — irY 
interaction. 
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FIG. 6: The total cross section of (a) K p — > K p, (b) K p — > ir + Y, (c) K p — > n S + (d) 
K~p — » 7T°S and (e) K~p — » 7r°A reactions in the relativistic model. The solid (dashed, dotted) 



curve shows the result using model (aA) ((aB), (fA)). Data are taken from Ref. 
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3, 



31 



32, 





A 


A(MeV) 


scattering length 


Relativistic 


0.90 


800 


0.175m- 1 


Nonrelativistic 


0.85 


800 


0.177m" 1 



TABLE V: Parameters and scattering length for the relativistic and nonrelativistic model of 1=1/2 
ttN interaction. 



The parameters of the model are A and a for the form factor and the strength parameter A. 
The obtained parameters are summarized in Table IVI1 The relativistic model can describe 
well the phase shifts up to 1.2 GeV as shown in Fig. [7J however, the non- relativistic model 
starts to deviate from the data at around 1.1 GeV. 

We used a Yamaguchi-type separable interaction for the nucleon-nucleon potential. To 
take into account the long range attractive interaction and the short range repulsion of the 
two-nucleon interaction, we used a two-term separable potential, 

< p\V BB \p>= C R g R (i?)g R (p) - C A g A (i?)g A (p). (27) 
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A 


A(MeV) 


a(fm) 2 


scattering length 


Relativistic 


2.7 


618 


0.50 


-0.095m" 1 


Nonrelativistic 


3.0 


628 


0.30 


-0.101m" 1 



TABLE VI: Parameters and scattering length for the relativistic and nonrelativistic model of 1=3/2 
ttN interaction. 
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FIG. 7: The phase shift of the irN scattering for Su (a), and S31 (b) partial waves. The solid 
curve shows the relativistic model and the dashed curve shows the nonrelativistic model. Data are 
taken from Ref. [.' 



Here Cr (Ca) is the coupling strength of the repulsive (attractive) potential. Quip) (9a{p)) 
is the form factor, whose form is given as gn{p) = A 2 R /(f^ + A^) (gA^P) = A- A /(jf + A^)), 
where A is a cutoff of the nucleon-nucleon potential. The adjustable parameters in our 
nucieon-nucleon potential are determined by fits to the data of the 1 Sq phase shifts 



The best-fit parameters are summarized in Table IVII The low energy phase shifts of the 1 Sq 
state is shown in Fig. [HI 





A R (MeV) 


A A (MeV) C R (MeVfm 3 ) 


C A (MeVfm 3 ) 


Relativistic 


1144 


333 


5.33 


5.61 


Nonrelativistic 


1215 


352 


5.05 


5.84 



TABLE VII: Our parameters of the relativistic and nonrelativistic model for NN scattering. 
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50 100 150 200 250 300 



T Lab (MeV) 

FIG. 8: Phase shifts of the NN scattering for 1 Sq state. Solid curve shows the relativistic model 
and dashed curve shows the nonrelativistic model. The phase shifts calculated from the model of 

n 

Ref. |36l ] are shown in triangles. 

IV. RESULTS AND DISCUSSION 

The dibaryon resonance with J 71 = 0~ , S = —1, 1=1/2 is studied using a formalism 
of the Faddeev equation as explained in section 2. We assume all the angular momentum 
to be in an s-wave state and the spin singlet state Sbb = for the two baryon states. We 
have included the dominant KNN, irUN and 7iAN Fock-space components, whose isospin 
wave functions are [K ® [NN] I=1 ] I=1 / 2 , [it® [SA r ]/=i/2,3/2]/=i/2 and [ir ® [AN] 1=1/2)1=1/2- An 
approximation within this model is that the weak YN interaction is not included. 

Let us start to examine the three-body resonance energy by taking into account only 
the KN interactions v^'^ N neglecting the ttYN Fock space. In this case, the bound 
state pole is expected to lie on the physical Riemann sheet below mx + 2rriN if the KN 
attraction is strong enough. Therefore it is not necessary to use the analytic continuation of 
the amplitude with the deformed contour discussed in section 2, so we simply use the integral 
over the momentum in the real axis. The results are shown in Fig. marked by a and a' 
for the the 'relativistic' and 'non-relativistic' models. Here we use the 'standard' parameters 
(aA) of the KN interaction with nonrelativistic and relativistic kinematics. The binding 
energies are about 18 MeV. The A^A" interaction included in r and Z is strong enough to 
bind the KNN system, where the 1 = KN interaction plays a dominant role. We then 
take into account the AW interaction. Then the binding energy increased furthermore to 
25.1 MeV (22.8 MeV) shown as b (b') for 'relativistic' ('nonrelativistic') model. Notice that 
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FIG. 9: The pole trajectories of the KNN — ttYN scattering amplitude for the J w = 0~ and 
1=1/2 state. The solid curve and the filled circles (dashed curve and filled triangles) show the 
results of the relativistic (nonrelativistic) model (aA). Here Wknn = rn-K + 2rriN- 

if we neglect the repulsive component of the NN interaction, we obtain a much more deeply 
bound state. 

In the next step, we gradually include the ttYN interactions, while the pion-exchange 
Z diagram is not yet included. To do this, we multiply by factor x the coupling constants 
Ca t p of the KN — ttY and 7rY — nY interactions as xC a ^. When the parameter is zero, 
x = 0, the ttY is disconnected from KN and when it takes the value 1, x = 1, we recover 
the full model. By varying the parameter x from to 1, we can follow the trajectory of 
the resonance pole from the bound state pole. Now the KNN bound state decays into 
the iiYN channel and the bound state pole moves into the unphysical sheet. Since the 
KNN bound state was found above the 7r£iV threshold, the resonance pole may be on the 
nYN unphysical and KNN physical Riemann sheet, which we have discussed in section 
2. The results of the pole trajectories are shown by the solid and dashed curves in Fig. [9] 
corresponding to the relativistic and the nonrelativistic models. Increasing the coupling to 
the 7rYN channel causes the width as well as the binding energy of the resonance increases. 
For larger binding energy Re(W po [ e — Wknn) < ~ 60 MeV the width starts to decrease 
because of the decreasing phase space for the decay into the 7r£iV state. The pole position 
is at — 82 — i29 MeV (—91 — 28i) for the relativistic (nonrelativistic) model shown as c(c'). It 
is noticed that the numerical method to follow the pole trajectories helps us to find whether 
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-100 -80 -60 -40 -20 



Re[W p0 | e -W KNN ] (MeV) 

FIG. 10: The pole trajectories of the three-body resonance (solid curve) and A(1405) (dashed 
curve) using the nonrelativistic model of (aA). 

we encounter singularities or not. As an example, the pole of the three-body resonance is 
shown by the solid line in Fig. [10] for < x < 1. The pole of the A(1405) is also shown by 
the dashed curve. The trajectory of the KNN resonance is similar to the the case A in Fig. 
IH and the integration contour does not intercept the singularity arising from the two-body 
resonance A(1405). 

Finally we include the 7r exchange mechanism in Z and ttN two-body scattering terms in 
r, which adds another mechanism for the decay of the KNN into tcYN and is important for 
the width of the three-body resonance. The final results of the KNN — nYN resonance poles 
are denoted by d and d' in Fig. [9j This mechanism increases the width of the three-body 
resonance by about 14 MeV, while the effect on the real part is small. The cancellation 
between the attractive 1 = 1/2 ttN interaction and the repulsive I = 3/2 irN interaction 
may lead to the small effects on the real part of the resonance energy. The effects of 7rAiV 
channel is small which increases binding energy and half- width at most by lMeV. The pole 
position of the three-body resonance is W = M — iY/2 = 2m^ + — 79.3 — £37.1 MeV 
(2rriN + rrix — 92.2 — i35.4 MeV) for relativistic (nonrelativistic) model shown as d and d'. 

The model dependence of our results on the three-body resonance is summarized in Tables 
IVIIII and IIXI The KNN — ttYN resonance pole is located on the KNN physical and ttYN 
unphysical sheet with the binding energy, B ~ 60 — 95 MeV, and the width, Y ~ 45 — 80 
MeV, using relativistic models. All of our models predict resonance energies above the ttY,N 
threshold. The relatively large model dependence of our results is due to the uncertainty 
in the models of I = KN — 7rE interaction. Comparing the results of model (A) with 



19 





Model (A) 


Model (B) 


(a) 


-79.3-«'37.1 


-79.3- i37.3 


(b) 


-93.3 - i27A 


-93.3-i27.6 


(c) 


-57.2-i38.6 


-56.9 - i38.6 


(d) 


-72.4-^31.7 


-72.2 - i31.9 


(e) 


-87.1-i40.8 


-87.1 - i41.0 


(f) 


-63.3 -»22. 2 


-63.2 - i22.3 



TABLE VIII: The pole energy (W po i e — mx — ^win) of the three-body resonance using relativistic 
models. The listed pole energies in MeV can be related to binding energy B and the width T as 
Wp i e - uik - 2m N = -B - iT/2. 





Model (A) 


Model (B) 


(a) 


-92.2 - i35.4 


-92.3-^35.6 


(b) 


-101.6 -»20.7 


-101.6 - z20.7 


(c) 


-72.7-i53.9 


-72.5 - «54.9 


(d) 


-83.0 - i33.3 


-83.0-i33.6 


(e) 


-98.1-i33.2 


-98.2-^33.3 


(f) 


-66.5 - i24.4 


-66.3 - «24.4 



TABLE IX: The pole energy of the three-body resonance. The same as Table IVIIII but for the 
nonrelativistic models. 

model (B), we can see the three-body pole position is almost independent of the parameters 
of the 1 = 1 KN — tiY interaction. By varying the real (imaginary) part of the fitted 
scattering length by ±0.1 fm, the binding energy of the three-body resonance is affected by 
~ ±14(8) MeV. Following another way to construct the model, the parameters of the model 
(f) are fitted to the pole energy of the A(1405). This model predicts the scattering length 
— 1.72 ± i0.44 fm. The energy of the three-body resonance is found to be B = 63 MeV with 
a rather small width, T = 44 MeV, compared with the models (a-e), which can already be 
seen in the small imaginary part of the scattering length in model (f). 

Let us briefly compare our results with those of the other theoretical studies of the K~pp 
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resonance, which use a nonrelativistic approach. Our resonance has a deeper binding energy 
and a similar width compared with those in Ref. [4]. However, it is not straightforward 
to compare with the pole energy of Ref. [4] because of the differences in the method to 
obtain the three-body resonance energy and on the model for the KN interaction. Their 
KN potential is stronger and has a short range than ours. Recently Shevchenko, Gal and 
Mares [3] studied K~pp system using the nonrelativistic coupled- channel Faddeev equation. 



Though the details of their method is not described in Ref. [19|, it seems their approach is 
quite similar to our present study. They employed a phenomenological KN potential model, 
and reported B ~ 55 — 70 MeV and T ~ 95 — 110 MeV. Their result is consistent with our 
results of the nonrelativistic model. Specially our result using the model (c) gives a quite 
similar resonance energy and width. 

In summary we have studied the existence and properties of a strange dibaryon resonance 
using the KNN — ixYN coupled channel Faddeev equation. By solving the three-body 
equation the energy dependence of the resonant KN amplitude is fully taken into account. 
The resonance pole has been investigated from the eigenvalue of the kernel with the analytic 
continuation of the scattering amplitude on the unphysical Riemann sheet. The model of 
the KN —ttY interaction is constructed from the leading order term of the chiral Lagrangian 
takes into account the relativistic kinematics. The KN interaction parameters are fitted to 
the scattering length given by Martin. We found a resonance pole at B ~ 79 MeV and 
r ~ 74 MeV in the relativistic model (aA). However, as the KN interaction is not very well 
constrained by the data, we studied a possible range of the resonance energies by considering 
different parameter sets of the KN — nY interaction. The binding energy and the full width 
can be in the range of B ~ 60 — 95 MeV and T ~ 45 — 80 MeV when computed in the 
relativistic model. In order to connect the resonance found in this work to the experimental 
signal, further theoretical studies on the production mechanism and further decay of the 
resonance especially to the A — p channel are necessary. 
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APPENDIX 



The spin-isospin recoupling coefficient of the particle exchange interaction Z is briefly 



explained. The coefficient given in Eq. (1.181) of Ref. [21J] can be simplified for s-wave 
states. The three-body state with total spin and isospin {S to tihot)i which couples with 
baryon 'isobar' with spin and isospin (S,I) and the spectator baryon B^SiM), is given as 

l[[ M 3(5 3 ,/3) ® Bj(SsJi)h,i) ® s Wo](W*.«) > • ( A - X ) 

Here baryon i,j represents particle 1 or 2, and the meson is always assigned as the third 
particle. The wave function of the three-body state, which couples with dibaryon 'isobar' 
and the spectator meson M3, is given as 

\[[Bi( Su h) ® B 2i s 2j2) } { s,i) ® M H s 3 ,i 3 )}(s tot ,itot) > • (A.2) 

Then Eq. (J7J) is extended to include spin-isospin degrees of freedom. The particle exchange 
interaction for the spectators l,m, the isobars /', / with spin-isospin (S',I') and (S,I) and 
the exchanged particle n can be expressed as follows 

^9f(s',i')(qi)g.f(s,i)(qm) 



Zl,f'{S',r),m,f{S,I)(PhPm,W) = Rl,f'(S',I'),m,f{S,I) J d(j>l ' Pr, 



W - Eiipi) - E m (p m ) - E n (pi + p m ) 

(A.3) 



where / represents isobar Y K , Y n , d and N*. 

Ri,f(S',i'),m,f{S,i) is given by the overlap of the initial and final spin-isospin wave functions. 
For the meson(M 3 ) exchange mechanism, Rij'(s',r)jj(s,i) is given as, 

Rij'(s',i'),jj(s,i) = < [[M 3 (s 3 ,i 3 ) ® £j(s,->/j)](sv') ® ^i(Si,/i)](s tot ,/ tot )l 

x |[[A f 3( < ff 8 ,/3) ® S W0](S,J) ® S j(5 J -,/ J )](S tot ,J tot ) > 



(_l)5+5'-5 3 -5 toW( ^ ; ^ 5tQt) 5 .. 5j lS ") v /( 2l S'+l)(2 1 S" + l) 

x (_i) i+i'-i*-i«*w{Ii, I 3 , I tot , If, I, I') v/(2/+l)(2/' + l). (A.4) 
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For the baryon(-Bj) exchange mechanism, Rij>(s',r),3,f(s,i) is given as, 

RiJ'(S',I'),3J(S,I) = < i[ M 3(S 3 ,I 3 ) ® Bj^.j.^S'jr) ® A(5 1 ,Ji)](5t«.t,Jtrt)l 
x |[[ 5 l(Si,/i) ® ^2(S 2) / 2 )](S,/) ® Af3(S 3 A)](S to t,/iat) > 

= (_ 1 )S 3 +s-s tot+ / 3+ /-/ toW((S , 3) ^ ^ 5 . . s ,/ j 5)^(25" + !) (25 + 1) 



x W{I 3 , Ij, hot, U\ I', I) v / (2/' + l)(2/ + l) 

x(Mi,i + M^-l)*^~ ff+Wi ~')- (A.5) 



When we anti-symmetrize the AGS equation, the last factor in the bracket in Eq. flA.5j) 



projects the anti-symmetric two nucleon states. This can be explicitly seen by comparing the 
exchange of nucleon 2 Z 2i Y K (s',i'),3,d(s,i) and nucleon 1 Z l>Ylc (S',i'),3,d(S,i) interactions. Using 



Eq. (1A.5I) . those interactions are related as 

Rl,Y K {S',I'),3,d{S,I) = (-i) S+I R2,Y K (S',I'),3,d{S,I), (A. 6) 

which leads to Eq. (JTT|) as 

XykXk = (1 - (-l) S+I )Z YK ,dTd,dXd,Y K + ■■■ ■ (A.7) 
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